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This paper introduces a general computational model for regenerators fed by a hot fluid stream on one 
side, whereas the other side is a fluid mixture, in which one of the components (refrigerant) undergoes 
a phase change. A simplified physical model, which combines fundamental and empirical correlations, 
and principles of classical thermodynamics, mass and heat transfer is developed and the resulting three- 
dimensional differential equations are discretized in space using a three-dimensional cell centered finite 
volume scheme. The proposed model is used to simulate numerically the behavior of the regenerator 
working under different operating and design conditions. Mesh refinements are conducted to ensure the 
convergence of the numerical results. The proposed methodology is shown to allow a coarse converged 
mesh for all simulations performed, therefore combining numerical accuracy with low computational 
time. The model was then used to determine optimal parameters for a specific geometric configuration of 
the regenerator, finding that the two-way maximized efficiency operating with optimal porosity and 
dimensionless mass flow rate of liquid ammonia, is Jjmax.max = 45.5% for the optimal pair 
(</>,m s ) opt = (0.5,0.01). As a result, the model is expected to be a useful tool for simulation, control, 
design, and optimization of regenerators for heat driven refrigerators. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Utility companies worldwide acknowledge that heating, 
ventilation, air conditioning and refrigeration (HVAC-R) systems 
are responsible for roughly 30% of total energy consumption, 
including the domestic, commercial and industrial sectors [1], 
therefore unquestionably with a major impact on energy demand. 
Researchers in many countries have been involved in developing 
refrigeration systems that deal with the drawbacks of conventional 
systems. 

Economic and environmental considerations brought a new 
point of view about refrigeration supplied by renewable heat 
sources. The international refrigeration industry has been investing 
considerable resources in that direction [2—13], For example, an 
absorption system or ejector cooling system driven by waste heat 
designed adequately could supply the requirement of oil platforms 
or refineries in regard to refrigeration and air conditioning. 

Absorption and ejector refrigeration systems can compete with 
compression technology in some applications. Also due to the use 
of cogeneration systems with turbine and gas engines in industrial 
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sectors, interest has increased in cooling machines which use 
residual heat [ 14]. Development of these machines requires high 
efficiency in the heat and mass transfer processes which take place 
in the absorber, which in turn is the most important component in 
these systems [15,16]. 

Constructal theory originated within the field of heat transfer 
[17,18] and has grown into a branch of thermodynamics that deals 
with engineering design. System configuration is obtained from 
minimization of resistance (to flow, to heat transfer). In this paper 
the regenerator is thermodynamically optimized as a result of the 
balance of competing factors along the lines of other previous 
studies [19]. 

A number of papers related to thermodynamic optimization of 
regenerator heat exchangers can be found in the literature [20—32]. 
Devices studied include rotary [20—24] or periodic-flow [25,26] 
regenerators, moving bed heat exchangers [27], and their applica¬ 
tions to gas turbines [28] and Stirling engines [29], as well as 
Stirling-cryogenic [30] and pulse tube refrigerators [31,32]. These 
optimization efforts involved the determination of the following 
optimal values and objective functions: optimum channel length 
and operating period for maximum second-law efficiency [20], and 
how the latter is affected by heat conduction in the regenerator 
matrix and unbalanced heat-capacity rates [21]; optimum opera¬ 
tional conditions for maximum thermal effectiveness of an air-to- 
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Nomenclature 

£ 

tolerance 




porosity 

A 

area, m 2 

V 

efficiency 

/Igl 

heat transfer area between matrix and hot fluid at the 

Pg 

density of metallic matrix, kg m~ 3 


cell “m”, m 2 

Pr 

density of fluid refrigerant, kg nT 3 

C 

specific heat, J kg -1 I< 1 

Pnh 3 

density of ammonia, kg itT 3 

C P 

specific heat at constant pressure of hot fluid, 

Ph 2 o 

density of water, kg m~ 3 


J kg^ 1 K 1 

P 

density of hot fluid, kg m~ 3 

Cv 

specific heat at constant volume of hot fluid, J kg -1 IC 1 

V 

specific volume, m 3 kg 1 

Cpr 

specific heat at constant pressure of internal fluid, 




J kg 1 IC 1 

Subscripts 

Cl 

ratio between the diameter of innermost tube to hot 

abs 

absorber 


fluid tube 

ar 

wire 

C 2 

ratio between the diameter of inner tube to hot fluid 

cond 

condenser 


tube 

cond.a 

conduction between metallic matrix of cell “m” and 

D 

diameter of external tube (hot fluid), m 


metallic matrix of cell “m - 1” 

d 

diameter, m 

cond.p 

conduction between metallic matrix of cell “m” and 

e 4 

thickness of innermost tube, m 


metallic matrix of cell “m + 1” 

h 

heat transfer coefficient by convection, W itT 2 K 1 

g 

metallic grid 

h t 

heat transfer coefficient by convection between cold 

gl 

interface between matrix and hot fluid 


fluid and internal tube, W m~ 2 IC 1 

gs 

transversal interface between two cells 

hf g 

enthalpy of vaporization of ammonia, J kg 1 

H 

high 

k 

thermal conductivity of metallic matrix, W m _1 K 1 

H 2 0 

water 

k s 

thermal conductivity of innermost tube, W m 1 K 1 

in 

inner tube or inlet 

L 

length of heat exchanger, m 

in,4 

initial state of liquid ammonia 

m 

mass, kg 

1 

liquid ammonia 

m 

mass flow rate, kg/s 

L 

low 

m 

dimensionless mass flow rate of hot fluid 

max 

maximum 

m s 

dimensionless mass flow rate of liquid ammonia 

nh 3 

ammonia 

TlceW 

number of cells 

NH 3 ,s 

liquid ammonia or vapor of ammonia completely 

Pin 

pressure of liquid ammonia that enters in the 


vaporized 


regenerator, bar 

NH 3 ,v 

vapor of ammonia vaporized accumulate 

Q 

heat transfer rate or enthalpy transferred, W 

opt 

optimum 

Qe nth3 

energy rate through the system 3 at cell “m", W 

pr 

internal fluid at constant pressure 

Qe nth4 

energy rate through the system 4 at cell “m", W 

r 

cold fluid 

Qin 

enthalpy transferred from cell “m” to cell “m - 1", W 

ref 

reference 

Q out 

enthalpy transferred from cell “m” to cell “m + 1", W 

S 

section interface of innermost tube 

Q 

heat transfer rate, W 

s 

liquid ammonia 

Qm 

heat transfer rate received by boiling ammonia, W 

s,0 

initial state of ammonia in the inner tube 

r 

fraction mass of water in strong solution 

s,in 

inlet state of ammonia in the inner tube 

T 

temperature, K 

sat 

saturated 

t 

dimensionless temperature 

T 

refrigerant inside of innermost tube 

t 

time, s 

T 

total 

u b 

ratio between global heat transfer coefficient between 

t 

internal tube 


system 3 and system 4 while phase changing to global 

tl 

interface between system 1 and system 3 


heat transfer coefficient 

tl4 

lateral interface of innermost tube 

U s 

global heat transfer coefficient between system 3 and 4 

V 

empty space in system 1 


when ammonia is liquid or vapor, Wm 2 IC 1 

V 

ammonia vapor 

U 34 

global heat transfer coefficient between system 3 and 

vr 

internal fluid at constant volume 


system 4, W irr 2 K 1 

w 

ammonia in the innermost tube 

U 

velocity of refrigerant, m s _1 

w,in 

inlet ammonia that does not vaporize 

V 

volume, m 3 

0 

initial state 

^ref 

length of reference, m 

0 

ambient 

X r 

quality of ammonia 

3 

interface between system 3 (non-vaporized ammonia) 

X™ 

quality of ammonia in the cell “m” 


and innermost tube 

v^cell 

A r 

quality of ammonia in the last cell 

4 

interface between innermost tube and liquid ammonia 

X 

dimensionless length 

00 

reference for temperature 

Greek symbols 

Superscripts 

Af 

amount of time that the fluid stays n the cell “m”, s 

m 

position of the cell (control volume) in the regenerator 

Ax 

length of cell “m", m 

ftcell 

last cell 
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air rotary regenerator [22], later extended to a multi-objective 
optimization (pressure drop and effectiveness); optimum length 
of the heat exchanger and core porosity for maximum heat transfer 
rate per unit of frontal surface area [24]; the geometry of the matrix 
of a periodic-flow heat exchanger with the unit cost of the exergy of 
the warm delivered air as the objective function [25]; geometry and 
operational parameters for minimum entropy generation [26]; 
optimal dimensions for maximum energy performance [27]; 
optimum temperature efficiency of the regenerator providing 
minimum heat transfer surface area [28]; design features to enable 
best Stirling engine efficiency [29]; and the effect of design 
parameters on the exergetic efficiency of the regenerator [30,32]. 
All these works rely, of course, on some sort of numerical modeling 
of the regenerator. In this literature review, no previous studies on 
the optimization of regenerators heat exchangers applied to 
absorption refrigeration systems were found. 

An experimental system was built in the laboratory using an 
absorption refrigerator and a car engine. The hot gas from the 
engine exhaust provides the waste heat used to drive the absorp¬ 
tion refrigerator. Fig. 1 shows a photograph of the prototype and 
Fig. 2 shows the corresponding schematic diagram illustrating the 
main system components and energy interactions. 

This paper introduces an unsteady mathematical model for 
a regenerative heat exchanger, which is part of heat driven refrig¬ 
eration systems that use waste heat (e.g., hot exhaust gases). The 
main objective is to investigate the existence of optimum operating 
and geometric parameters for maximum system efficiency. 


m,T n 


Absorption 

Refrigerator 



cond 


Condenser 


Absorber 


Evaporator 



m,T in 


Regenerator 
Heat Exchanger 


— 


§* 


T, 


Fig. 2. Schematic diagram with the main components of the absorption refrigerator 
equipped with a regenerator heat exchanger. 


2. Mathematical model 

The regenerator consists of a heat exchanger in which a metallic 
matrix is mounted inside the hot fluid passage to provide energy 
storage capabilities. The metallic matrix is particularly important 
when the flow rate of hot fluid is not enough to keep the system 
working at the desired levels of performance. In these situations, 
the amount of energy (heat) saved by the metallic matrix will 



sustain the amount of heat transferred to the cold fluid despite of 
the low hot fluid mass flow rate. 

The amount of metallic matrix inside the hot fluid tube is rep¬ 
resented in the model by the porosity, <p = V v /W- 

2.1. Cold fluid without phase change 

As a preliminary step, it is studied how the amount of metallic 
matrix affects the performance of the regenerator when there is no 
phase change of the cold fluid. 

Fig. 3 shows the regenerator divided into small control volumes 
named “cells”, in which L = n ce nAx. 

Each cell of the regenerator is divided in three different systems. 
System 1: metallic matrix and tube; system 2: hot fluid, and system 
3: cold fluid. Fig. 4 shows the thermal interactions among the 
systems in a generic cell “m” of the regenerator. 

2.1.1. System 1; metallic matrix and tube 

Applying the first law of thermodynamics to system 1, the 
resulting equations are 


Qg + Qcond.p ' Qcond.a ' Qt 


dr- 

= dt 


(i) 


Cold flow 


Hot 

flow 



Ax 


Fig. 1. Prototype of absorption refrigerator equipped with a regenerator heat 

exchanger. Fig. 3. Heat exchanger divided into control volumes. 
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Adiabatic 



nr» m+1 
g 


where q g = h g A$(T m - T g m ); q CO nd,p = -M gs ((T g m - T g m+1 )/Ax), 
qcond,a = -fcA gs ((T g m - T™ _1 )/Ax), and q t = MfP(7'r m - T™). 

Porosity is an important parameter in the model because it 
affects the thermal response of the regenerator. It is possible to 
express the areas in Equation (1) as functions of porosity, as 


follows: 


A t ] = ttdxl 

(2) 

„ 4(1 - 0)V T 

Agl ” dar 

(3) 

A gs L = (1 - $)Vj 

(4) 

where Vj = (jt/ 4)((D 2 - d 2 )L). The regenerator is insulated, there¬ 
fore when m = 1, q C ond,a = 0- an d when m = n ce h, qcond.p = 0. 

In this study, the heat transfer coefficients by convection, h g and 
fit are taken as constants, but it is possible to account for their 
dependence on the Reynolds (Re) and Prandt (Pr) numbers using 
empirical correlations [33]. 

The total amount of mass of metallic matrix at cell “m" is 
expressed as a function of total volume of metal in the regenerator, 
as follows: 

— Pg^metal 

where VWtai = A gs L. 

(5) 

2.2.2. System 2: hot fluid flow 

The energy balance for system 2 is 


m dl m 

Qin - <7out -Qg = m m Cv-g^ 

(6) 

where q in = rfiCpT 111-1 and q ou t = rhc p T m . 

The mass of hot fluid in the mth cell, 
a function of porosity as follow 

m m , is expressed as 

m m = pcpVj 

(7) 


When m = 1 (first cell), T 111 1 = Tj n . 


2.13. System 3: cold fluid flow 

The energy balance for system 3 is 

-q t + m r c pr (^T™ +1 — T r m ) = m ™ c vr^b (8) 

where = p r (Ttd 2 /4)L 

Equations (1), (6) and (8) comprise a system of three differential 
equations to be integrated in time for each cell. Considering that 
n ce n is the number of cells, 3 n ce ii is the total number of ordinary 
differential equations to be integrated in time for the entire 
regenerator. 

Next, the possibility of cold fluid (refrigerant) phase change in 
the regenerator is included in the analysis. 

2.2. Cold fluid with phase change 

Fig. 5 shows the thermal interactions among the systems that 
occur in any cell “m” of the regenerator when the cold fluid expe¬ 
riences a phase change. In the end of the regenerator, the refrig¬ 
erant (ammonia) has to be totally vaporized in order to go forward 
in the cycle; if the vaporization of ammonia is not complete, the 
reminder of liquid ammonia shall return to the beginning of the 
regenerator and pass through it again. Because of these features, 
one more flow is inserted in the regenerator, as shown in Fig. 5. 

Hence, the model includes hereafter four different systems, i.e.: 
system 1: metallic matrix and tube; system 2: hot fluid; system 3: 
cold fluid in the inner tube (non-vaporized ammonia) and system 
4: cold fluid in the innermost tube (liquid and vapor ammonia 
mixture). 

Ammonia and impurities (water) flow through the innermost 
tube. Phase change occurs as a result of heat received from the hot 
fluid flow. Depending on geometric and operating conditions some 
amount of liquid ammonia might not vaporize. The liquid ammonia 
that did not vaporize, returns through the inner tube (system 3) to 
a reservoir located before the regenerator in order to re-enter the 
regenerator. 

In the model, the water present in ammonia is assumed to be 
just a small fraction of the mixture (e.g., ejector refrigeration 
systems [8]), so that the refrigerant (ammonia) is treated a pure 
substance. If considered as a binary mixture, the temperature 
during phase change would no longer be constant, and the 
temperature variation should follow specific binary mixtures phase 
diagrams [34,35], 
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Fig. 5. Heat and mass transfer at cell “m” when there is phase changing of cold fluid (ammonia). 


2.2.1. System 1: metallic matrix and tube 

The equations developed for system 1 with no phase change of 
cold fluid can be used to model system 1 with phase change. This is 
possible because the mass and heat transfer that occur in system 1 
do not depend on the ammonia phase change. 

2.2.2. System 2: hot fluid flow 

The energy balance for system 2 is 

rhCp^T m_1 - T m J - q s = m m c v -gp (9) 

when m = 1 (first cell), T" 1-1 = Tj n . 


during phase change, a ratio between the global heat transfer 
coefficient when there are two phases in the system (vapor and 
liquid) to when just one phase is present (vapor or liquid), Ub was 
defined as follows: 


Uu = 


U 34 

U s 


(H) 


The mass of liquid ammonia (if present) in system 3 at cell “m” 
can be expressed as a function of the quality of ammonia as follows: 


‘NH, 


= (1 -< cdl )(l -r)r 


( 12 ) 


where r is the mass fraction of water in the strong solution defined 
as r = rii H2 o/m T and m T = ih NH3 , s + m H2 o- 

The total mass of refrigerant in cell “m” in the system 4 is 


2.2.3. System 3: cold fluid flow in the inner tube (non-vaporized 
ammonia) 




\xA s 

v 


(13) 


The energy balance for system 3 is 

—<Jt + <Js + <Jenth 3 = ( m NH 3 c NH3,l + m H 2 O c H 2 o) pjp 
where 

<7s = U 34 A™{T™-T™)-,q en th3 

= (mNH3 C NH 3J + rilH 20 CH 2 o) ' (C.in “ C) 


where v = (1 - r)u NH3 + ru H2 o- 

The mass of water in ammonia at cell “m” is 

( 10 ) o = rmf (14) 


2.2.4. System 4: cold fluid flow in the innermost tube (liquid 
ammonia) 


and 


The energy balance for system 4 is 


itfdts + ^L 
tW ^cel ' 

U 34 can assume two different values depending on whether 
there is phase change or not. When there is no phase change: 


/ N dT m 

<7enth4-<7s = (mNH 3 ,s c NH 3 ,i + m[? 2 oCH 2 oJ -gp '= 0, V) (15) 

While liquid ammonia is vaporizing the temperature is kept 
constant and equal to the saturation temperature. Hence, in this 
case: 


U34 = U 5 


1 

± e4 ± 

h 3 + k s + h 4 


While ammonia is vaporizing the global heat transfer coefficient 
assumes numerical values greater then U s . In order to express U 34 


djm 

^ = 0 
df 


(16) 


Equation (16) is applied from the first cell “m” where the 
vaporization begins until the phase change is complete (or the last 
cell of the regenerator is reached). 
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In order to choose between Equations (15) and (IS) to model the 
behavior of system 4, it is necessary to find out where the phase 
change starts and finishes or where there is no phase change. This 
can be done checking if Tf” > T sa t and if the mass of vapor of 
ammonia generated according to the flow is less or equal than the 
total mass of ammonia at cell “m”, i.e., m}J H3 v < (1 - r)mY- If these 
two conditions are true, then the refrigerant at cell “m” is changing 
phase, otherwise the refrigerant is liquid (in the beginning of the 
regenerator) or vapor in the end of the regenerator. 

The temperature TY is calculated solving the Equation (15) or 
assuming that TY = Tsat during phase change. 

As long as TY > T sat , v will remain constant and equals to 
zero. The mass of refrigerant that evaporates in each cell “m” is 


Am 


m 

NH 3 ,v 


-QsAt 

frfg(Tsat) 


(17) 


where At — Ax/u and u — m T v/A s . 

The mass of refrigerant vapor accumulated in cell “m” is 


m 


m 

NH 3 ,v 


m 


m 

NH 3 ,v 


+ Am 


m 

NH 3 ,v 


(18) 


The next step is to calculate the quality of the regenerator in cell 
'm", xJT 


x 


m 

r 


111 


m 

NH 3; v 


(1 - r)mf 


(19) 


When m[J H3 v > (1 - r)my (or x™ > 1), the phase change is 
complete and Equation (16) does not express the behavior of 
system 4 anymore, then Equation (15) shall be used. 

The mass of liquid ammonia or ammonia totally vaporized in 
cell “m” and the mass flow rate of ammonia in cell “m” are 
respectively m™ H3 s = (1 - r)mf and m NH3 , s = (1 - r)m T . 

The energy rate through system 4 in the cell “m” is 


<7enth4 = (m[MH 3 ,sCNH 3 ,i + ™H 2 oCH 2 o) ' (T^ - T s m ) (i = l,v) 

( 20 ) 

It is important to emphasize that in this model, the refrigerant is 
considered a pure substance because the mass fraction of water in 
the refrigerant is very small (r = 0.01). This value of mass fraction 
allows for the assumption that the temperature during phase 
change remains constant and equal to the ammonia saturation 
temperature. 

The energy balance for system 4 is expressed by Equation (15) 
when there is one phase in the system (liquid or vapor) and it is 
expressed by Equation (16) when there is phase change. 


2.3. Efficiency 


The efficiency (tj) of the regenerator is defined as the ratio 
between the actual heat transfer rate to ammonia vapor and the 
maximum possible heat transfer rate. 


v _ Q nh 3 ,v 


(21) 

Qmax 


where 



9max = rilCp(r in -r ini 4) 


(22) 

<Jnh 3 ,v = (1 - r)x? ce "m T \ 

[ c NH 3 ,l(^sat — Tj n4 ) + hfg(Tsat) 

+ CNH 3 ,v(T^'-r sat )} 

(23) 



x 

Fig. 6. Temperature distribution in the regenerator for ti> = 0.95. 


The regenerator heat exchanger is modeled according to two 
different conditions. In the case where there is no phase change, 
the regenerator is represented in each cell by a system of three 
ordinary differential equations, Equations (1), (6) and (8); where 
TJ 1 , T™ 0 and Tfy are the initial conditions. In the case where the 
phase change is present, the regenerator is represented in each cell 
by a system of four ordinary differential equations, Equations (1), 
(9), (10), (15) or (16); where TJ 1 , T™ 0 , T™ 0 and are the 
initial conditions. For both cases, the systems were solved using 
a Runge—Kutta 4th/5th order method [36]. 



Fig. 7. Temperature distribution in the regenerator for 0 = 0.5. 
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x 

Fig. 8. Temperature distribution in the regenerator for </> = 0.2. 



t [s] 


Fig. 10. Temperature as a function of time in the first ceil when the hot fluid mass flow 
rate is modeled by a periodic function (0 = 0.95). 


The temperature of each system is obtained until the steady 
state is reached and after that the efficiency is calculated. The 
steady state is reached when the norm of the solution vector 
becomes independent of time according to a pre-specified toler¬ 
ance (e = 10 -4 ), as follows: 



(25) 


I r(t + At) - r'(t)H 
l|T(t)|| 


(24) 


and the number of cells in the mesh is taken as n ce iij when the norm 
of the steady state solution vector does not vary between a less 
refined (j) and a more refined mesh (j + 1) 


3. Results and discussion 

3.2. Model without cold fluid phase change 

The initial conditions for this configuration are To = 373.15 K, 
T r ,o = 293.15 K, m = 0.1 kg s _1 and L = 1 m. The numerical 



0 

Fig. 9. Time to reach steady operation according to different porosities (</>). 



Fig. 11. Temperature as a function of time in the first cell when the hot fluid mass flow 
rate is modeled by a periodic function (0 = 0.5). 
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Table 1 

Physical properties and geometric data used in 
the numerical solution of the systems. 

c g = 8961kg- 1 ic 1 

c p = 1000 J kg- 1 K 1 

c pr = c vr = 4180Jkg- 1 K- 1 

c v — 713 J kg -1 K- 1 

d ar = 0.002 m 

e 4 = 0.001 m 

h 3 = h 4 = lOOWm- 2 K- 1 

k = 14.9 W m- 1 1<-’ 

p = 1165 kg m- 3 

p % = 2707 kg rrr 3 


simulation of the regenerator is performed by solving Equations (1), 
(6) and (8). The scales = 298.15 K, L re f = 1 m, and 
rh re f = 1 kg s _1 are used to nondimensionalize the temperatures, 
length of the regenerator and mass flow rates, respectively, so that 
7j = Tj/Too (i = 1,2,3.4), x = x/L ref and m = m/m re{ . 

Figs. 6—8 show the numerically obtained dimensionless 
temperatures distribution in the regenerator heat exchanger for 
c/> = 0.95, 0.5, and 0.2, respectively. It is seen that for ^ = 0.2 the 
curves are closer to each other than for <j> = 0.95 and 0.5. This 
suggests that the heat transfer is more efficient when the porosity 



Fig. 12. quality distribution in the regenerator for r/t = 0.5. 


is small (more metallic mass), since the more metallic mass the 
more heat transfer area results between system 1 and system 3, as 
it is physically expected. 

The time to reach steady state operation as a function of porosity 
is shown in Fig. 9. The time to reach steady operation increases as 
the amount of metallic mass increases (porosity decreases) due to 
increasing thermal inertia, showing a trade off that has to be 
balanced in actual regenerator heat exchanger design. 

When the curves in Fig. 8 are too close to each other to a point in 
which they can be represented by a single line, it is reasonable to 
assume that there is no significant amount of heat being transferred 
in that specific region of the heat exchanger, therefore it could be 
shortened. So, for the specific configuration shown in Fig. 8, the 
heat exchanger could be 0.4 m long instead of 1 m long. This feature 
is very important when space is constrained, as for example, in 
automobile air conditioning systems. 

The next step is to study how the porosity affects the efficiency 
when the mass flow rate of hot fluid varies (due to some distur¬ 
bance). In order to simulate this behavior a periodic function was 
selected to define the mass flow rate. 

m = ^cos(t) + ^p (26) 

For this analysis, the simulation was conducted until t = 800 s. 
The results are shown in Figs. 10 and 11, for <f> = 0.95 and 0.5, 
respectively. The goal is to determine the behavior of the regener¬ 
ator when a disturbance occurs in the hot fluid mass flow rate. As it 
was concluded previously in the text, when the porosity decreases 
the temperature of the systems get close to each other due to 
a better thermal contact, however there is an additional and posi¬ 
tive effect. When the porosity is low, the effect of the disturbance is 
better “absorbed” by the metallic mass than with a higher porosity. 
Therefore the regenerator thermal response is more robust as 
porosity decreases, i.e., the system becomes less subject to 
perturbations such as fluctuations in heat supply, for example. 

3.2. Model with cold fluid phase change 

The initial conditions for this part of the analysis are the same as 
those used when there is no cold fluid phase change, but the 



Fig. 13. Efficiency as a function of dimensionless hot fluid mass flow rate for <j> = 0.5. 
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Fig. 14. Efficiency as a function of dimensionless hot fluid mass flow rate for different 
porosities. 

modeling equations are different. The numerical simulation of 
the regenerator is performed by solving Equations (1), (9), (10), and 
(15) or (16). 

The geometric ratios cj = d ts /D and C 2 = dj„/D are defined. The 
results for the physical properties and geometric data listed in 
Table 1 are plotted in Figs. 12 and 13, for 0 = 0.5, c\ = 0.3 and 
C 2 = 0.5. In Fig. 12, the variation of ammonia quality, x r , with 
respect to regenerator dimensionless length, x, is represented as 
a function of different liquid ammonia dimensionless mass flow 
rates, rh s . It is noted that, for rh s > 0.009 the ammonia flow does 
not vaporize completely, as it is physically expected since the heat 
input rate is the same as mass flow rate increases. Fig. 13 shows 
that the efficiency displays a maximum for rh s ,opt~0.01, i.e., the 
configuration with higher efficiency is the one that totally vapor¬ 
izes the ammonia flow and uses the entire heat exchanger length 
for that, delivering ammonia as saturated vapor at the outlet, i.e., 
neither as superheated vapor nor as a two-phase mixture (liquid 
and vapor). 

The main objective of the study is to find possible optimal 
configurations leading to the most efficient regenerator. Therefore, 
besides the dimensionless liquid ammonia mass flow rate, in 
Fig. 14, the effect of porosity on efficiency is investigated. The effi¬ 
ciency of the regenerator is depicted as a function of m s for different 
values of porosity, when Ci = 0.3 and C 2 = 0.5. According to Fig. 14 
the now two-way maximized efficiency is 45.5% for the optimal 
pair (0, rhs) opt = (0.5,0.01). 

The existence of the optima found with respect to m s in 
Fig. 14 is explained by the analysis of two extremes. For a fixed 
heat input rate, when the mass flow rate of ammonia is very 
large (approaches infinity), liquid ammonia mass flow rate to be 
heated is too high, thus the difference between inlet tempera¬ 
ture and outlet temperature is small, hence the actual heat 
transfer rate to ammonia vapor, qNH 3 ,v. is also small. When the 



<t> 

Fig. 15. Optimal hot fluid mass flow rate and efficiency as functions of porosity. 



Fig. 16. Efficiency of the regenerator for different values of (c 2 = 0.5). 


mass flow rate of ammonia approaches zero, q NH3 v is also small, 
because it is proportional to m s . Analogously, when 0 —> 0, the 
empty area available for the hot fluid to flow decreases so that 
the velocity of hot fluid increases. This effect leads to inlet and 
outlet temperatures with very similar values making the effi¬ 
ciency to decrease. When 0 —> 1, the heat transfer area between 
the hot fluid and metallic matrix decreases, so that the efficiency 
also decreases. 

Fig. 15 summarizes the two-way optimization procedure 
showing that the one-way maximized efficiency, i) max . with respect 
to m s , depicts a second maximum with respect to porosity, 0. 
Furthermore, m s ~ 0.01 for all porosities, i.e., the optimal operating 
condition is “robust” with respect to porosity variation, which is an 
important feature for regenerator operation. 

3.3. Parametric analysis 

Fig. 16 shows the results obtained when the porosity is kept 
constant (0 = 0.5), C 2 = 0.5 and the values of ci are varied. It is 
observed that there is no optimal ci, i.e., the efficiency increases as 
C\ increases. This feature occurs because in the model, the pressure 
drop was neglected in the pipes of the regenerator. In the optimi¬ 
zation of the geometric configuration of the regenerator the pres¬ 
sure drop is therefore an important feature to be considered. 

4. Conclusion 

In this paper, a mathematical model for a regenerative heat 
exchanger was developed based on the volume element method¬ 
ology to capture expected physical trends. Next, transient simula¬ 
tions were conducted and a thermodynamically optimized 
regenerator heat exchanger design was proposed. 

This study demonstrated that the insertion of a metallic matrix 
(represented by a porosity parameter in the model) helps the 
system to become more stable to possible disturbances in the mass 
flow rate supply of hot fluid. 

The model was then used to determine optimal parameters for 
a specific geometric configuration of the regenerator, finding that 
the two-way maximum efficiency operating with optimal porosity 
and dimensionless mass flow rate of liquid ammonia, is 
ifmax.max = 45.5% for the optimal pair (0,m s ) opt = (0.5,0.01). The 
main conclusion is that it is expected that, after experimental 
validation and model adjustment to ensure the accuracy of results, 
the developed mathematical model could be used as an efficient 
tool for the thermal design, control and optimization of regenerator 
heat exchangers with similar characteristics to the heat exchanger 
studied in this paper. 
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